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Abstract 

Background: Gastrointestinal bacteria play a central role in the health of animals. The bacteria that individuals 
acquire as they age may therefore have profound consequences for their future fitness. However, changes in 
microbial community structure with host age remain poorly understood. We characterised the cloacal bacteria 
assemblages of chicks and adults in a natural population of black-legged kittiwakes (Rissa tridactyla), using 
molecular methods. 

Results: We show that the kittiwake cloaca hosts a diverse assemblage of bacteria. A greater number of total 
bacterial OTUs (operational taxonomic units) were identified in chicks than adults, and chicks appeared to host a 
greater number of OTUs that were only isolated from single individuals. In contrast, the number of bacteria 
identified per individual was higher in adults than chicks, while older chicks hosted more OTUs than younger 
chicks. Finally, chicks and adults shared only seven OTUs, resulting in pronounced differences in microbial 
assemblages. This result is surprising given that adults regurgitate food to chicks and share the same nesting 
environment. 

Conclusions: Our findings suggest that chick gastrointestinal tracts are colonised by many transient species and 
that bacterial assemblages gradually transition to a more stable adult state. Phenotypic differences between chicks 
and adults may lead to these strong differences in bacterial communities. These data provide the framework for 
future studies targeting the causes and consequences of variation in bacterial assemblages in wild birds. 

Keywords: Age-differences, Automated ribosomal intergenic spacer analysis, Bacteria, Black-legged kittiwakes, 
Cloaca, Gastrointestinal tract 



Background 

In animals, bacterial communities occur both externally 
(e.g. skin, feather, scales) and internally (e.g. gastrointes- 
tinal and reproductive tracts) of their hosts [1]. Although 
many gastrointestinal bacteria are host commensal, 
others are adaptive by being central in fundamental 
physiological processes [2-4], or have pathogenic effects 
[5]. Given the potentially profound impact of the gastro- 
intestinal microbiota on host condition, the study of 
bacterial assemblages within animals is of high eco- 
logical and evolutionary importance [6]. 
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One of the central questions in the study of microbial 
communities in animals is how they are acquired and how 
the structure of communities varies with age. Young 
animals either passively or actively acquire microbiota 
from the environment after birth or hatching [6-8]. The 
gaxstrointestinal tract of young individuals is typically 
characterised by a high number of transient microbial spe- 
cies before eventually transitioning to a relatively stable 
state as adults [e.g. 9]. However, it remains unclear how mi- 
crobial communities change throughout the lifetime of 
their hosts and how these fluctuations affect host health, 
behaviour and fitness [6]. Temporal variation in microbial 
composition may have important consequences for factors 
such as nutrition e.g. [8,10], mating preferences e.g. [11] 
and susceptibility to diseases e.g. [12]. It is therefore neces- 
sary to characterise bacterial communities at different life 
history stages of animal hosts. Scant information exists, 
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however, and much past research has focussed on humans 
or other mammal species [7,8]. 

Birds comprise a useful model for studying the causes 
and consequences of variation in bacterial assemblages. 
For example, the contrasting life histories of birds and 
mammals may result in differences in the mechanisms 
of microbial colonisation. Mammals, for example, ac- 
quire important maternal microbes during the birth 
process [9], whereas birds are more likely to acquire mi- 
crobes after hatching from other sources, such as the 
nesting environment or food [7]. Second, many bird spe- 
cies regurgitate food to their young, thus permitting a 
mode of vertical microbial transmission that most mam- 
mals lack. Last, birds possess a cloaca, which serves the 
dual functions of excretion and gamete transfer. Gastro- 
intestinal bacteria may thus be transferred during copu- 
lation [13], providing the potential for pathogens [14,15] 
or beneficial strains to be sexually transmitted [16]. 

The few avian studies on the effect of age on gastro- 
intestinal assemblages have mostly been conducted on 
domestic birds e.g. [17-20]. However, the artificial rear- 
ing conditions in which captive animals are raised make 
it difficult to extrapolate knowledge of bacterial assem- 
blages in captive animals to their wild counterparts 
[21,22]. Moreover, the few studies in wild birds typically 
examine only a few target bacterial species and do not 
use community-level approaches to test for host age dif- 
ferences in bacterial assemblages e.g. [23-26], but see 
[10]. Very little is therefore known about how the nor- 
mal gastrointestinal bacterial communities change be- 
tween young and adult wild birds. This information 
forms an important basis for understanding bacterial ac- 
quisition and transmission mechanisms, the prevalence 
of enteropathogens and the consequences of particular 
communities of bacteria for host health. 

Here, we characterise the bacterial assemblages of both 
adults and chicks in a wild population of black-legged kit- 
tiwakes {Rissa tridactyla), a pelagic gull with a wide distri- 
bution. We focus our analyses on cloacal assemblages 
because they are known to approximate assemblages 
deeper within the gastrointestinal tract [27] while requir- 
ing relatively non-invasive sampling methods. We have 
previously used a molecular technique [ARISA - Auto- 
mated Ribosomal Intergenic Spacer Analysis: [28] to 
explore variation in cloacal bacteria assemblages in kitti- 
wakes [13]. Although ARISA allows the rapid assessment 
of bacterial assemblage diversity, it does not permit spe- 
cies identification. We have therefore further developed a 
comprehensive ribosomal RNA (rRNA) clone library to 
identify actual bacterial species. Combined with ARISA, 
this allows the rapid assessment of actual species in bac- 
terial assemblages, thus providing a framework for study- 
ing the causes and consequences of microbial assemblage 
variation in kittiwakes. We then use a community-based 



approach to test whether age predicts bacterial assemblage 
composition in kittiwakes and which bacterial lineages dif- 
fer between chicks and adults. 

Methods 

Cloacal bacterial assemblages were sampled from 22 adult 
and 21 black-legged kittiwake chicks captured at nests in 
an abandoned U.S. Air Force radar tower on Middleton Is- 
land (59° 26' N, 146° 20' W) in the Gulf of Alaska [13]. As 
an aim of our study was to taxonomically classify ARISA 
peaks (see below), it was necessary to obtain a comprehen- 
sive sample of the most common bacterial species present 
in chicks and adults. We therefore only sampled unrelated 
individuals (to minimise the overlap in bacterial assem- 
blages between two individuals due to, for example, a 
shared nesting environment). In addition, chicks of a wide 
range of ages were sampled (range = 5-30 days, mean = 
16.0 ± 8.3 SD days) as bacterial assemblages were 
hypothesised to change rapidly with age. Sampling oc- 
curred between May and July 2009. Cloacal sampling and 
DNA extraction protocols are outlined in detail in White 
et al. 2010 [13]. We also regularly collected control samples 
during field work by pipetting 1 ml of saline solution into a 
sterile vial from the saline stock, which ensured that all 
bacteria identified were of cloacal origin and not due to 
contamination of the stock solution or of reagents used for 
DNA extraction. 

Automated ribosomal intergenic spacer analysis 

We first characterised cloacal bacterial assemblages by 
performing PCR-ARISA, which exploits the extreme in- 
terspecies variability in the length of the intergenic spa- 
cer (IGS) lying between the conserved 16S and 23S 
rRNA genes in the bacterial ribosomal operon. ARISA 
involves the amplification of DNA extracted from the 
bacterial assemblage of interest using a fluorescently- 
labelled primer and subsequent high-resolution electro- 
phoresis in an automated system. Assemblages are 
therefore characterised by a series of electrophoretic 
peaks that vary according to the length of the amplified 
IGS fragment of each bacterial species. Given that unre- 
lated species may share the same IGS fragment length, 
the number of species represented by each ARISA peak 
remains uncertain. This method has been optimised for 
the characterisation of cloacal bacterial assemblages in kit- 
tiwakes [13]. Briefly, the 16S-23S rRNA intergenic spacer 
was amplified using the FAM (6-carboxyfluorescein)-la- 
beled primer S-D-Bact-1522-b-S-20 (5'-[6FAM] TGCG 
GCTGGATCCCCTCCTT-3') and the unlabelled L- 
D-Bact-132-a- A-18 (5'-CCGGGTTTCCCCATTCGG-3'). 
In addition to amplifying the highly variable IGS region, 
the primers also amplify a 131 bp fragment of the 23S 
rRNA gene, which we used for subsequent phylogenetic 
analyses (see below). 
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Fragment amplification was performed using PCR in 
20 ul reactions containing 0.2 mM dNTPs, 1.5 mM MgCl 2 , 
0.25 uM of each primer, and 1 unit of GoTaq DNA poly- 
merase with lx buffer (Promega). Due to the coextraction 
of host and environmental DNA, it was not possible to ac- 
curately measure the concentration of the extracted DNA. 
Therefore, a standard of 4 ul of eluted DNA extract was 
used per reaction. PCR was conducted using the following 
program: 94°C for 10 minutes, followed by 30 cycles of 
94°C for 1 minute, 55°C for 1 minute and 72°C for 1 mi- 
nute, and a final elongation step at 72°C for 10 minutes. 
The successful amplification of the samples was confirmed 
on a 2% agarose gel and all PCR products were then di- 
luted by a ratio of 1 in 30. One microlitre of the diluted 
product was then added to 8.7 ul Hidi Formamide (Applied 
Biosystems) and 0.3 ul of in-house designed ROX size 
standard (size standard range = 105 - 1007 bp). The mix- 
ture was then run on an ABI PRISM 3130x1 automated se- 
quencer (Applied Biosystems). Traces were viewed on 
GeneMapper software 4.0 (Applied Biosystems). 

IGS-23S rRNA library construction and analysis 

We created an IGS-23S rRNA clone library of our sam- 
ples in order to taxonomically classify ARISA peaks. We 
defined three hierarchical levels of units in our analyses: 
ARISA peaks, clones and OTUs (operational taxonomic 
units). As stated above, an ARISA peak is defined as one 
electrophoretic peak in an ARISA profile. Cloned frag- 
ments in our IGS-23S rRNA library were referred to as 
clones. During each cloning event we typically isolated 
multiple copies of the same sequence. Each unique clone 
(i.e. those with a unique sequence) was assumed to cor- 
respond to one unique ARISA peak. Finally, when mul- 
tiple unique clones shared the same 23S rRNA sequence 
and differed from each other only within the IGS, we as- 
sumed that they belonged to the same OTU. This ap- 
proach in OTU classification is conservative but justified 
given that many bacterial species have multiple operons 
in their genomes that result in multiple ARISA peaks 
per OTU [29]. 

In developing our rRNA library, IGS-23S fragments of 
our bacterial samples were amplified as outlined above, 
except that both primers were unlabelled. Amplified prod- 
ucts were then cloned using a TOPO TA cloning kit 
(Invitrogen) as per the manufacturer's protocol. We then 
picked between 8 and 96 transformed colonies for each 
cloacal sample and amplified the inserts using colony PCR 
in 25 ul reactions containing 2.5 mM MgCl 2 , 0.1 mM 
dNTPs, 0.2 uM M13 forward and reverse primer, and 2.5 
units of Firepol polymerase and lx buffer (Solis Bio Dyne). 
Reactions were run at 94°C for 5 minutes, followed by 35 
cycles of 94°C for 30 s, 55°C for 30 s and 72°C for 30 s, 
and a final extension at 72°C for 10 min. The successful 
amplification of the samples was confirmed on a 2% 



agarose gel and primers and excess dNTPs were removed 
from the amplified products by digestion with exonuclease- 
shrimp and alkaline phosphatase (Fermentas Life Sciences). 
Cloned PCR products were amplified in both directions 
using Big Dye chemistry (Applied Biosystems) and se- 
quenced on an ABI PRISM 3130x1 automated sequencer 
(Applied Biosystems). Sequence editing and alignment was 
conducted using CLC DNA Workbench 6.1 (CLC bio). 
Rarefaction analyses were conducted using Analytic Rar- 
efaction 1.3 (available at http://strata.uga.edu/software) to 
ascertain the comprehensiveness of our cloning effort. 

Comparison of ARISA peak size and clone fragment lengths 

One aim of our study was to identify the bacterial OTUs 
present in each cloacal sample of chicks and adults 
based entirely on the ARISA output. We therefore com- 
pared the length of cloned sequences with the peaks 
present in the ARISA profiles. We assigned each clone 
for each individual to an ARISA peak, taking into ac- 
count the presence or absence of the ARISA peak and 
clone sequence in the other cloned samples for that age 
class. Any peak to which we did not assign a sequence 
was assumed to be an artefact and was therefore elimi- 
nated from subsequent analyses. 

Bacterial taxonomy 

A Basic Local Alignment Search Tool search BLAST: 
[30] was implemented on the 23S rRNA sequence of 
cloned fragments to screen for published sequences that 
closely match our cloned sequences. The maximum 
BLAST score (typically > 200) and Expect-value (typic- 
ally < 1 x 10" 50 ) were used as indicators of the strength 
of sequence similarity. 

As BLAST provides information on similarity between 
two sequences but not on relatedness among sequences, 
we employed phylogenetic analyses to infer the taxo- 
nomic classification of the clones. We included all 
unique clones as well as two to four similar sequences 
for each clone from known species identified by BLAST. 
Both maximum-likelihood and Bayesian methods were 
used to ascertain the phylogenetic position of each un- 
identified species relative to known bacterial species. 
Given the immense diversity among bacteria, we 
conducted analyses separately for each phylum identified 
(Firmicutes, Actinobacteria and Proteobacteria). It was 
necessary to further separate the Proteobacteria into 
classes (Alphaproteobacteria, Betaproteobacteria and 
Gammaproteobacteria). Maximum-likelihood trees were 
constructed in Treefinder [31] using the nucleotide sub- 
stitution model suggested by the Bayesian information 
criterion. Bootstrap support was estimated after 1000 
replicates. Bayesian phylogenies were constructed in 
MrBayes 3.1.2 [32] using the same substitution model as 
suggested by Treefinder. Posterior probabilities were 
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estimated after 10 million generations, sampling every 
1000 generations. 

The most derived phylogenetic position for each clone 
was then estimated based on a combination of the BLAST 
result and the phylogenetic trees. For example, when both 
BLAST and the phylogeny agreed that a particular clone 
belonged to a specific genus, we classified the clone to that 
genus. However, if there was disagreement between 
BLAST and the phylogenetic tree in the phylogenetic pos- 
ition of the clone, we classified the clone to a higher taxo- 
nomic level (e.g. class or family), with an emphasis on the 
position of the clone within the tree. Full sequences for 
each OTU have been deposited in GenBank under Acces- 
sion numbers KC462595 - KC462734. 

Bacterial assemblage comparisons 

To characterise differences in the bacterial assemblages 
between chicks and adults, we implemented Unifrac [33] 
to compare overall bacterial assemblage composition 
using phylogenetic information. All assemblage analyses 
were conducted on the ARISA data and without abun- 
dance weights. Using the cloning data may produce mis- 
leading results due to stochastic variation in the identity 
of sequences cloned for each individual. 

A 23S rRNA phylogenetic tree containing all sequences 
found in chicks and adults was used for the analyses. We 
then implemented Unifrac to test 1) whether bacterial as- 
semblages cluster by age (i.e. adults and chicks) using both 
the environmental clustering analysis (setting 'number of 
sequences to keep' to five and 'number of permutations' to 
100) and principal coordinates analysis, 2) whether any 
clustering is statistically significant using both the P test 
(which does not use phylogenetic information; 100 permu- 
tations) and the Unifrac test (which tests whether more 
evolutionary history of bacterial assemblages is shared 
within age classes than between classes; 100 permutations) 
and 3) which bacterial lineages are responsible for differ- 
ences in assemblages between the age classes, via the 
lineage-specific analysis (setting 'minimum descendants' 
to four). 

A greater quantity of bacterial DNA could be extracted 
from adults than chicks (data not shown), probably due to 
the larger volume of saline solution that could be used to 
flush adult cloacae. Cloning of bacterial DNA fragments 
therefore tended to be more successful for adults than 
chicks and twice as many clones were sequenced per indi- 
vidual for adults compared to chicks (number of clones se- 
quenced: adults - 47.0±28.0 clones, chicks - 20.2±6.7 
clones, Mann Whitney U = 434.0, P<0.001). This age- 
related discrepancy could lead to biases in downstream 
analyses in, for example, our estimates of age differences 
in the number of species hosted per individual. This, how- 
ever, does not appear to be the case. First, the proportion 
of ARISA peaks that were discarded because they could 



not be assigned OTUs did not differ between adults and 
chicks (proportion of peaks discarded per individual: 
adults - 0.63±0.13, chicks - 0.66±0.14, F 1>41 =0.534, 
P=0.469), although the absolute number discarded was 
higher in adults (number of peaks discarded per individ- 
ual: adults - 12.7±5.9 peaks, chicks - 7.6±1.6 peaks, Mann 
Whitney U = 101.0, P=0.001) probably because higher 
quantities of bacterial DNA in adults produced more noise 
in the ARISA profiles. Second, cloning effort did not differ 
between chicks of different ages, both when correlating 
chick age with number of clones sequenced per individual 
(r = 0.273, F lil9 =1.528, P=0.231) or in a univariate test 
grouping chicks into young (5-10 days old) and old indi- 
viduals (20-30 days old; number of clones sequenced per 
individual - young chicks: 18.5±4.5 clones, old chicks: 
22.4±8.3 clones, F U9 =1.966, P=0.177). Last, our Unifrac 
analyses are unaffected by these potential biases as we 
kept a minimum of five sequences for the clustering ana- 
lyses and ran 100 permutations on the data. All non- 
genetic statistical analyses were conducting using SPSS 
17.0 (SPSS, Chicago, Illinois, USA). All data are presented 
as mean ± standard deviation. 

Results 

Comparison of ARISA peak size and clone fragment 
lengths 

A mean of 33.9 ± 24.4 clones per individual were se- 
quenced (range = 8 - 120 clones per individual). Clones 
ranged in size from 304 bp to 1090 bp (Table 1). In total, 
we identified 142 unique clones, of which 99 were suc- 
cessfully matched to ARISA peaks. In the remaining 43 
cases, we were unable to match peaks to clones because 
the clone was rare, resulting in low confidence of assign- 
ment. Due to the relatively high proportion of OTUs for 
which we could not assign ARISA peaks, the analyses 
presented below (e.g. differences in bacterial assemblages 
between adults and chicks) apply only to the most com- 
mon OTUs found in chicks and adults. 

The 142 unique clones corresponded to 76 OTUs. This 
discrepancy was due to indels causing variation in clone 
lengths for several OTUs. For example, although some 
OTUs were always consistent in clone length (e.g. OTU 2, 
Table 1), others were more variable (e.g. OTU 5 varied be- 
tween 342 and 345 bp in length). In most cases, one 
ARISA peak corresponded to a single OTU. However, in 
16 cases, one ARISA peak was assigned to two distantly- 
related OTUs. For example, the 348 bp peak corresponded 
to species within the genera Corynebacterium and 
Athromitus, while the 544 bp peak corresponded to species 
of Rothia and Staphylococcus. The mean number of ARISA 
peaks per OTU was 1.5 ± 1.5 (range 1-9). Of the 76 OTUs 
identified, 64 were of probable cloacal origin. The 
remaining 12 were also isolated from control samples and 
so were assumed to be contaminants (Additional file 1). 
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Table 1 Identity of OTUs isolated from cloacae of black-legged kittiwakes 



OTU 


ARISA 
peak 


Clone 


Difference 


r>i a c t j- _ i - 

BLAST match 


BLAST 
score 


Most derived 
phylogenetic position 


Phylum 


Found 


1 


304 


308 


+4 


Peptoniphilus asaccharolyticus 


152 


Genus: Peptoniphilus 


Firmicutes 


Ad 


! 


305 


309 


+4 


Peptoniphilus asaccharolyticus 


152 


Genus: Peptoniphilus 


Firmicutes 


Ad 


! 


396 


399 


+3 


Peptoniphilus asaccharolyticus 


152 


Genus: Peptoniphilus 


Firmicutes 


Ad 


! 


397 


400 


+3 


Peptoniphilus asaccharolyticus 


152 


Genus: Peptoniphilus 


Firmicutes 


Ad 


, 


398 


401 


+3 


Peptoniphilus asaccharolyticus 


152 


Genus: Peptoniphilus 


Firmicutes 


Ad 


2 


313 


315 


+2 


Enterococcus pseudoavium 


228 


Genus: Enterococcus 


Firmicutes 


Ch 


3 


331 


331 


0 


Enterococcus faecal is 


211 


Order: Lactobacillales 


Firmicutes 


Ad, Ch 


3 


432 


433 


+1 


Enterococcus faecal is 


211 


Order: Lactobacillales 


Firmicutes 


Ad, Ch 


4 


340 


341 


+1 


Enterococcus malodoratus 


239 


Genus: Enterococcus 


Firmicutes 


Ch 


4 


341 


342 


+1 


Enterococcus malodoratus 


239 


Genus: Enterococcus 


Firmicutes 


Ch 


4 


433 


434 


+1 


Enterococcus malodoratus 


239 


Genus: Enterococcus 


Firmicutes 


Ch 


5 


345 


342 


-3 


Corynebacterium propinquum 


239 


Genus: Corynebacterium 


Actinobacteria 


Ad 


5 


346 


343 


-3 


Corynebacterium propinquum 


239 


Genus: Corynebacterium 


Actinobacteria 


Ad 


5 


347 


344 


-3 


Corynebacterium propinquum 


239 


Genus: Corynebacterium 


Actinobacteria 


Ad 


5 


348 


345 


-3 


Corynebacterium propinquum 


239 


Genus: Corynebacterium 


Actinobacteria 


Ad 


6 




342 




Enterococcus malodoratus 


206 


Order: Lactobacillales 


Firmicutes 


Ch 


/ 


348 


349 


+1 


Arthromitus sp. 


233 


Genus: Arthromitus 


Firmicutes 


Ad, Ch 


/ 


418 


419 


+1 


Arthromitus sp. 


233 


Genus: Arthromitus 


Firmicutes 


Ad, Ch 


8 




349 




Lactobacillus salivarius 


206 


Genus: Lactobacillus 


Firmicutes 


Ch 


9 


349 


350 


+1 


Lactobacillus animalis 


213 


Genus: Lactobacillus 


Firmicutes 


Ad, Ch 


9 


440 


443 


+3 


Lactobacillus animalis 


213 


Genus: Lactobacillus 


Firmicutes 


Ad, Ch 


9 


441 


444 


+3 


Lactobacillus animalis 


213 


Genus: Lactobacillus 


Firmicutes 


Ad, Ch 


9 


519 


520 


+1 


Lactobacillus animalis 


213 


Genus: Lactobacillus 


Firmicutes 


Ad, Ch 


9 


520 


521 


+1 


Lactobacillus animalis 


213 


Genus: Lactobacillus 


Firmicutes 


Ad, Ch 


10 


352 


353 


+1 


Aerococcus urinae 


222 


Order: Lactobacillales 


Firmicutes 


Ad 


1 1 


357 


354 


-3 


Corynebacterium aurimucosum 


196 


Genus: Corynebacterium 


Actinobacteria 


Ad 


12 




357 




Atopobium parvulum 


204 


Order: Clostridiales 


Firmicutes 


Ad 


13 




357 




Carnobacterium mobile 


237 


Genus: Carnobacterium 


Firmicutes 


Ch 


13 




463 




Carnobacterium mobile 


237 


Genus: Carnobacterium 


Firmicutes 


Ch 


14 




357 




Lactobacillus crispatus 


237 


Genus: Lactobacillus 


Firmicutes 


Ad 


14 


603 


607 


+4 


Lactobacillus crispatus 


237 


Genus: Lactobacillus 


Firmicutes 


Ad 


15 




361 




Bacillus megaterium 


215 


Genus: Bacillus 


Firmicutes 


Ch 


16 


365 


361 


-4 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


16 


366 


362 


-A 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


16 


367 


363 


-4 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


16 


368 


364 


-4 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


16 


369 


365 


-4 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


16 


371 


367 


-4 


Mycobacterium sp. 


161 


Suborder: Corynebacterineae 


Actinobacteria 


Ad 


17 




363 




Oceanobacillus iheyensis 


187 


Class: Bacilli 


Firmicutes 


Ch 


18 


380 


375 


-5 


Arcanobacterium abortisuis 


1/--1 


Family: Actinomycetaceae 


Actinobacteria 


Ad 


19 




37/ 




Clostridium difficile 


217 


Genus: Clostridium 


Firmicutes 


Ch 


20 


376 


378 


+2 


Enterococcus faecal is 


239 


Genus: Enterococcus 


Firmicutes 


Ch 


20 


478 


480 


+2 


Enterococcus faecal is 


239 


Genus: Enterococcus 


Firmicutes 


Ch 


21 


392 


388 


-4 


Corynebacterium aurimucosum 


182 


Genus: Corynebacterium 


Actinobacteria 


Ad 
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Table 1 Identity of OTUs isolated from cloacae of black-legged kittiwakes (Continued) 



22 


388 


390 


+2 


Pediococcus inopinatus 


174 


Order: Lactobacillales 


Firmicutes 


Ch 


23 


392 


395 


+3 


Clostridium piliforme 


219 


Genus: Clostridium 


Firmicutes 


Ch 


23 


393 


396 


+3 


Clostridium piliforme 


219 


Genus: Clostridium 


Firmicutes 


Ch 


24 




398 




Acidimicrobineae sp. 


180 


Phylum: Actinobacteria 


Actinobacteria 


Ad 


25 




400 




Clostridium piliforme 


209 


Genus: Clostridium 


Firmicutes 


Ad 


25 


753 


756 


+3 


Clostridium piliforme 


209 


Genus: Clostridium 


Firmicutes 


Ad 


26 


398 


401 


+3 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


399 


402 


+3 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


406 


40/ 


+1 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


430 


431 


+1 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


431 


432 


+1 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


489 


492 


+2 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


520 


520 


0 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 


521 


521 


0 


Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


26 




621 




Staphylococcus vitulinus 


243 


Genus: Staphylococcus 


Firmicutes 


Ch 


27 


421 


422 


+1 


Bacillus sp. 


213 


Genus: Bacillus 


Firmicutes 


Ch 


27 


422 


423 


+1 


Bacillus sp. 


213 


Genus: Bacillus 


Firmicutes 


Ch 


27 


570 


568 


+2 


Bacillus sp. 


213 


Genus: Bacillus 


Firmicutes 


Ch 


27 


589 


590 


+1 


Bacillus sp. 


213 


Genus: Bacillus 


Firmicutes 


Ch 


28 




434 




Streptococcus parauberis 


231 


Genus: Streptococcus 


Firmicutes 


Ch 


29 




438 




Acidimicrobiaceae sp. 


189 


Phylum: Actinobacteria 


Actinobacteria 


Ch 


30 


437 


439 


+2 


Staphylococcus saprophyticus 


239 


Genus: Staphylococcus 


Firmicutes 


Ch 


30 


438 


440 


+2 


Staphylococcus saprophyticus 


239 


Genus: Staphylococcus 


Firmicutes 


Ch 


31 


446 


441 


-5 


Corynebacterium propinquum 


182 


Genus: Corynebacterium 


Actinobacteria 


Ad 


32 


443 


445 


+3 


Bacillus globisporus 


206 


Genus: Bacillus 


Firmicutes 


Ch 


33 




446 




Bacillus cereus 


213 


Genus: Bacillus 


Firmicutes 


Ch 


34 


452 


448 


-4 


Corynebacterium propinquum 


172 


Genus: Corynebacterium 


Actinobacteria 


Ad 


34 


458 


454 


-4 


Corynebacterium propinquum 


172 


Genus: Corynebacterium 


Actinobacteria 


Ad 


34 


459 


455 


-4 


Corynebacterium propinquum 


172 


Genus: Corynebacterium 


Actinobacteria 


Ad 


34 


460 


456 


-4 


Corynebacterium propinquum 


172 


Genus: Corynebacterium 


Actinobacteria 


Ad 


34 


461 


45/ 


-A 


Corynebacterium propinquum 


172 


Genus: Corynebacterium 


Actinobacteria 


Ad 


35 


453 


454 


+1 


Corynebacterium amycolatum 


200 


Genus: Corynebacterium 


Actinobacteria 


Ad 


35 


454 


455 


+1 


Corynebacterium amycolatum 


200 


Genus: Corynebacterium 


Actinobacteria 


Ad 


36 


480 


479 


-1 


Rhodococcus erythropolis 


207 


Genus: Rhodococcus 


Actinobacteria 


Ad 


36 


490 


490 


0 


Rhodococcus erythropolis 


207 


Genus: Rhodococcus 


Actinobacteria 


Ad 


36 


493 


492 


-1 


Rhodococcus erythropolis 


207 


Genus: Rhodococcus 


Actinobacteria 


Ad 


37 


502 


505 


+3 


Escherichia coli 


213 


Genus: Escherichia 


Proteobacteria 


Ad, Ch 


3/ 


504 


50/ 


+3 


Escherichia coli 


213 


Genus: Escherichia 


Proteo bacteria 


Ad, Ch 


37 


506 


509 


+3 


Escherichia coli 


213 


Genus: Escherichia 


Proteobacteria 


Ad, Ch 


37 


596 


599 


+3 


Escherichia coli 


213 


Genus: Escherichia 


Proteobacteria 


Ad, Ch 


37 


598 


601 


+3 


Escherichia coli 


213 


Genus: Escherichia 


Proteobacteria 


Ad, Ch 


38 




508 




Bacillus globisporus 


198 


Genus: Bacillus 


Firmicutes 


Ch 


39 




509 




Nocardia seriolae 


1/4 


Order: Actinomycetales 


Actinobacteria 


Ch 


39 




514 




Nocardia seriolae 


1/4 


Order: Actinomycetales 


Actinobacteria 


Ch 


40 


520 


522 


+2 


Corynebacterineae sp. 


207 


Order: Actinomycetales 


Actinobacteria 


Ad 
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Table 1 Identity of OTUs isolated from cloacae of black-legged kittiwakes (Continued) 



A 1 


527 




_1 


IVIILI UUULLC.I lUl ! 1 LC.bLuLt.UI I / 


1 °>8 


OfcrllUi. IVIIL! UUULLC.I lU! ! 1 


ML LI I lUUdLLfcrl Id 


Ch 


Al 
4Z 


S.A1 
J4Z 


^3A 


T 
— Z 


Rothio mucilaginoso 


1 7 A 
1 /4 


Genus: Rothio 


Acti no bacteria 


AH Th 


42 


544 


538 


_T 


PiUttHU 1 1 lULIIUyll lUbU 


1 74 


Con ic" Qr^thi'tfi 


ALU I lUUdLLfcrl Id 


Ad Ch 


43 


544 


545 


+ 1 


J>LU{JI lyiULULLUj UUic.Ub 


20 7 


UtrllUb. J>LU{J! lyiULULLUi 


I 1 1 1 I IILULfcrb 


Ad 


A A 




^zia 

J^+O 




Nocordia s&noiQB 


1 

I JU 


Order! Actinomycetales 


Acti no bacteria 




4j 




Jj4 




Loctobociilus mucosae 


zuy 


Genus: Loctobociilus 


Firmicutes 


AH 


/1A 
40 


c; 7t 

J/Z 


J// 


I c 
tj 


Clostridium tyfobutyficum 


Of"! 

yu 


Order: Clostridiales 


Firmicutes 


AH 
AO. 


A 1 
4/ 




_>uU 




Streptococcus pyogenes 


1 OQ 

i yo 


Genus: Streptococcus 


Firmicutes 


ah r\~\ 
au, 


A / 

4/ 


/U/ 


/l 1 
/ 1 z 


, c 

TJ 


Streptococcus pyogenes 


1 OQ 

i yo 


Genus: Streptococcus 


Firmicutes 


ah r[~, 
Au, Ln 


A / 

4/ 


/ 1 j 


71 Q 
/ I 0 


+ J 


Streptococcus pyogenes 


1 OQ 

i yo 


Genus: Streptococcus 


Firmicutes 


ah r[~, 
AO, Ln 


47 


/ 1 6 


721 


+5 


Streptococcus pyogenes 


1 98 


Genus: Streptococcus 


Firmicutes 


Au, Ln 


47 


785 


788 


+3 


Streptococcus pyogenes 


1 98 


Genus: Streptococcus 


Firmicutes 


Au, Ln 


40 


COT 

jyz 


COI 

jy i 


— i 


Micfococcus luteus 


T1 3 
Z I Z> 


Genus: Micrococcus 


Actinobacteria 


Ln 


4y 


ai 1 

O I I 


ouy 


— Z 


Corynebacterium efficiens 


1 AO 

i oy 


Genus: Corynebacterium 


Actinobacteria 


Ln 


_)U 


£1 1 


ai ^i 

D 14 


_T_3 


Streptococcus dysgaloctiae 


T1 3 
Z I J 


Genus: Streptococcus 


Firmicutes 


AH 
AU 


c; 1 

J I 




A1 C 
0 I J 




Jonesio denitrificons 


1 QT 

I Oz 


Suborder: Micrococcineae 


Actinobacteria 


rh 
Ln 


j i 




OZU 




Jonesio denitrificons 


1 RT 


Suborder: Micrococcineae 


Actinobacteria 


rh 


j i 




OZ I 




Jonesio denitrificons 


1 QT 

I OZ 


Suborder: Micrococcineae 


Actinobacteria 


rh 
Ln 


j l 




ATT 




Jonesio denitrificons 


1 QT 
I OZ 


Suborder: Micrococcineae 


Actinobacteria 


rh 

LI l 


^i 


ATT 
0ZZ 


atc 
DZj 


+J 


Sinobacteraceae sp. 


T 1 1 
Z 1 1 


Order: Xanthomonada es 


Proteo bacteria 


AH 
AU 


jj 




ATA 

uzu 




Agrococcus sp. 


in 7 
zU/ 


Family: Microbacteriaceae 


Actinobacteria 


rh 
Ln 


J4 




0j4 




Ciavibocter michigonensis 


1 Qt; 

I O J 


Family: Microbacteriaceae 


Actinobacteria 


rh 

LI I 


cc 




a 7n 
o/u 




Yersinio enterocolitica 


Tit; 

Z 1 J 


Genus: Yersinio 


Proteo bacteria 


AH 
AU 


c;a 
JO 




A /I 
O/ 1 




Herbospiriilum sp. 


T 1 / 
Z 1 / 


Genus: Herbospiriilum 


Proteo bacteria 


AH 
AU 


j/ 




O/ 




Arthrobocter oriioitensis 


1 OR 

i yo 


Family: Micrococcaceae 


Actinobacteria 


rh 


Jo 




AQT 

oyz 




Azoorcus sp. 


1 cm 
i y j 


Family: Rhodocyclaceae 


Proteo bacteria 


rh 
^_i i 


59 




713 




Agromyces sp. 


185 


Family: Microbacteriaceae 


Actinobacteria 


Ch 


60 




723 




Micrococcus luteus 


176 


Family: Micrococcaceae 


Actinobacteria 


Ch 


61 




740 




Arthrobocter oriioitensis 


182 


Family: Micrococcaceae 


Actinobacteria 


Ch 


62 


794 


796 


+2 


Phenylobacterium zucineum 


165 


Family: Caulobacteraceae 


Proteobacteria 


Ad 


63 


973 


974 


+1 


Bradyrhizobium elkonii 


209 


Genus: Bradyrhizobium 


Proteobacteria 


Ad, Ch 


63 




1047 




Bradyrhizobium elkonii 


209 


Genus: Bradyrhizobium 


Proteobacteria 


Ad, Ch 


63 




1048 




Bradyrhizobium elkonii 


209 


Genus: Bradyrhizobium 


Proteobacteria 


Ad, Ch 


64 




1090 




Bradyrhizobiaceae sp. 


198 


Genus: Bradyrhizobium 


Proteobacteria 


Ch 



For each OTU, the following variables are displayed: OTU identification number, ARISA peak size, clone fragment length, the difference between the latter two 
variables (i.e. clone length minus ARISA peak size), the species identity of the most similar sequence in GenBank as determined by the BLAST algorithm, the 
maximum BLAST score, the most derived phylogenetic position of the OTU and the phylum to which the OTU belongs. In some cases, we could not determine 
which ARISA peak represented a particular clone (indicated by a blank space in the "ARISA peak" column). 



A virtually perfect correlation existed between the size 
of ARISA peaks and clone length (r = 1.000, Fi, 95 = 299,083, 
P < 0.001), although some discrepancies still existed (Table 1). 
Clone lengths differed from their corresponding ARISA peak 
sizes by a mean of 1.0 ± 3.6 bp (range -6 to +15 bp). 

Bacterial taxonomy 

The resolution of OTU identification depended on both the 
strength of the BLAST match and the ability to resolve their 



respective phylogenetic positions (Additional file 2). The 
mean maximum BLAST score for each OTU was 199 ± 25 
(range = 90 - 243). In 64.5% of cases we assigned an OTU 
to a specific genus, to a family in 14.5% of cases, to a sub- 
order in 2.6% of cases, to an order in 13.2% of cases, to a 
class in 2.6% of cases and to a phylum in 2.6% of cases. The 
most common phylum represented by the cloacal assem- 
blages was Firmicutes (50.0%), followed by Actinobacteria 
(37.5%) and Proteobacteria (12.5% consisting of 4.7% 
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Alphaproteobacteria, 3.1% Betaproteobacteria and 4.7% 
Gammaproteobacteria). The most common genera found 
were Corynebacterium (10.9%) and Lactobacillus (6.6%). 
The majority of OTUs identified in the control samples 
belonged to the class Alphaproteobacteria (58.3%). In ad- 
dition, we also identified bacteria from the classes Beta- 
proteobacteria (16.7%) and Gammaproteobacteria (8.3%) 
and from the phylum Firmicutes (16.7%). 

Bacterial assemblages in adults and chicks 

From our 23S rRNA cloning library, we identified 31 
OTUs in adults and 40 OTUs in chicks, despite sequen- 
cing less than a quarter the number of chick to adult 
clones. Rarefaction analyses revealed that we had sam- 
pled the majority of common OTUs found in adults, but 
that many rare OTUs remained unidentified in chicks 
(Figure 1). Only seven OTUs were shared between the 
two age groups. Sixty three percent of OTUs identified 
in chicks were isolated from only single individuals 
while, in adults, 48% of OTUs were isolated from single 
individuals. From the ARISA data, chick age was posi- 
tively correlated with number of OTUs hosted per indi- 
vidual (r = 0.452, F 1;19 = 4.889, P = 0.039; Figure 2). In 
addition, adults hosted a higher number of OTUs per indi- 
vidual than chicks (mean number of OTUs per individual 
kittiwake: chicks = 5.8 ± 3.0 OTUs, adults = 9.2 ± 2.3 
OTUs, F li43 = 17.556, P < 0.001). 

Analyses in Unifrac revealed strong clustering of the 
sequences by age (Figure 3). In addition, principal coor- 
dinates analysis separated chicks and adults along the 
first axis explaining 45.6% of the variation (PCI scores - 
chicks: 0.264 ± 0.070, adults: -0.252 ± 0.091, F lj41 = 
430.5, P < 0.001; Figure 4), but not along the second 
axis, which explained 12.0% of variation, or the third 
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Figure 1 Rarefaction curves of OTU richness and cloning effort 
of cloacal bacteria of adult and chicks. Adults are represented by 
the dashed line and chicks by the solid line. Dotted lines represent 
95% confidence limits. 



axis, explaining 7.9% of variation (PC2 scores - chicks: - 
0.019 ± 0.195, adults: 0.018 ± 0.099, Z = -0.875, P = 
0.382; PC2 scores - chicks: -0.019 ± 0.195, adults: 
0.018 ± 0.099, Z = -0.462, P = 0.664). 

The clustering by age was statistically significant (P test 
in UniFrac: P = 0.007), although the clusters did not repre- 
sent evolutionary divergence of bacteria between the age 
classes (Unifrac significance: P = 0.104). Lineage-specific 
analysis suggested that the age differences in bacterial as- 
semblages were due to the presence of several OTUs in 
only one age class. For example, various OTUs of the 
genus Corynebacterium were abundant in adults (six 
OTUs identified, many of which were present in most in- 
dividuals), but virtually absent in chicks (one OTU identi- 
fied in only one individual; P < 0.001). Similarly, certain 
OTUs were only found in adults (e.g. OTU 1: 
Peptoniphilus sp., P < 0.001; OTU 18: Actinomycetaceae 
sp., P < 0.001; OTU 46: Clostridiales sp., P < 0.001), while 
others were unique to or more abundant in chicks (e.g. 
OTU 27 and 32; Bacillus sp., P < 0.001; OTU 37: 
Escherichia sp., P < 0.001). 



Discussion 

We identified 64 bacterial OTUs in the cloacae of adult 
and chick black-legged kittiwakes, a majority of which 
were identified as Firmicutes and Actinobacteria, while a 
smaller proportion belonged to the Proteobacteria. A 
greater number of total OTUs were identified in chick clo- 
acae than in adults. However, the number of OTUs hosted 
per individual increased with age. Older chicks hosted 
more OTUs per individual than younger chicks and adults 
hosted more OTUs per individual than chicks. Surpris- 
ingly very little overlap existed in the bacterial assemblages 
between chicks and adults. Only seven of 64 OTUs were 
shared between the two age groups (e.g. OTU 3: order 
Lactobacillales, OTU 9: genus: Lactobacillus and OTU 37: 




10 20 
Chick age (days) 

Figure 2 Age-related changes in the number of bacterial OTUs 
hosted by black-legged kittiwake chicks. 
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Figure 3 Unifrac clustering analysis of bacterial assemblages in 
the cloacae of adult and chick black-legged kittiwakes. Adults 
are represented by the dark grey box and chicks by the light grey box. 
Each leaf represents the bacterial assemblage of one individual 
kittiwake. Nodes with numbers were recovered more than 95% of the 
time during 100 permutations. Note that the analysis removed several 
chick individuals that harboured fewer than five OTUs (see Methods). 



genus Escherichia), while the vast majority were found ex- 
clusively in only chicks or adults. For example, Corynebac- 
terium was the most common genus of bacteria identified 
in adults, but it was virtually absent from chicks. This pro- 
nounced difference between adults and chicks resulted in 
strong statistical clustering of bacterial assemblages 
according to host age. 

ARISA has been widely used to provide important in- 
sights in a wide range of fields within microbiology e.g. 
[34-38]. Coupled with a clone library for species identifi- 
cation, it represents a highly relevant research tool to 
allow the rapid and inexpensive characterisation of en- 
vironmental bacterial assemblages. Despite the advan- 
tages of ARISA, this technique has some limitations 
(which are not all necessarily restricted to ARISA). For 
example, biases inherent during DNA extraction and 
PCR are also known to affect the apparent composition 
of bacterial assemblages [39,40]. ARISA can potentially 
underestimate species richness, as eight percent of bac- 
terial species are known not to have the 23S and 16S 
rRNA organised in an operon (i.e. they have no 
intergenic spacer region) or have very large IGS lengths 
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PC1 

Figure 4 First two principal coordinates extracted from a 
principal coordinate analysis of kittiwake cloacal bacterial 
assemblages. Chicks are represented by closed circles and adults 
by open circles. The first principal component explained 45.6% of 
the variation, while the second explained 12.0%. 



that cannot be detected by ARISA [29]. In addition, di- 
vergent bacteria may share the same IGS length and 
therefore be associated with the same ARISA peak [this 
study, [29]. In contrast, some species have several op- 
erons in their genomes resulting in multiple ARISA 
peaks for single species and potentially leading to over- 
estimates in diversity [29]. In our dataset, we conserva- 
tively assumed that any two clones with the same 23S 
sequence derived from the same OTU. We consequently 
identified many OTUs with multiple operons (1-9 op- 
erons were identified per OTU). We also only included 
OTUs in our analyses for which we could confidently as- 
sign ARISA peaks, resulting in several genuine OTUs be- 
ing excluded. This conservative approach means that our 
estimates of OTU richness are likely to represent a mini- 
mum. However, these biases apply equally for both chicks 
and adults and we were still able to identify a relatively 
large number of OTUs suitable for community-level 
comparisons. 

The establishment of bacterial communities in the 
gastrointestinal tract of young animals is characterised 
by a high turnover of many transient species and large 
changes in community structure over short periods of 
time. For example, Gonzalez-Braojos et al. 2012 [23] 
found that, in faecal sacs of nestling pied flycatchers 
{Ficedula hypoleuca), Enterobacteriaceae loads decreased 
when the nestlings aged from 7 to 13 days, while Entero- 
cocci loads concurrently increased. Age-related changes 
in the composition of other important gastrointestinal 
microbes, such as fungi, are also known to occur e.g. 
[25] . The rapid changes in bacterial community structure 
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in young animals may arise due to a number of reasons in- 
cluding resource competition between bacterial species, 
shifts in host diet or age-related variation in the chemical 
and physiological state of the gastrointestinal tract 
[8,10,23]. For example, in humans, the early colonisation 
of the gut by facultative anaerobes (e.g. Enterococci and 
Enterobacteria) reduces gut oxygen levels which allows 
anaerobic bacteria to become established [8]. Eventually, 
gastrointestinal bacterial communities are known to tran- 
sition to a stable adult state [7]. Adult bacterial communi- 
ties may differ from those of young individuals due to the 
more developed immune system of adults [41-43], the low 
mobility of young animals resulting in a restricted envir- 
onment from which to obtain bacteria or contrasting 
chemical and anatomical cloacal environments that are 
differentially hospitable or hostile to various bacteria [7,8]. 

Our data support these findings. For example, the fact 
that the chick rarefaction curve failed to plateau and that 
more OTUs were unique to individual chicks than adults 
suggests that chicks host more rare, and presumably 
transient, bacterial species than adults. Second, our find- 
ing that adults host a greater number of OTUs per indi- 
vidual than chicks, and that the number of OTUs hosted 
by chicks increases with age, is in accordance with previ- 
ous studies that have shown that species richness in bac- 
terial assemblages increases as animals reach adulthood 
e.g. [19,20,26]. Last, we identified substantial variation in 
bacterial assemblage composition between chicks and 
adults, highlighting the dynamic nature of bacterial com- 
munities within the gut. It is, unfortunately, not possible 
to deduce the fitness consequences on hosts of age- 
related changes in bacterial microbiota from our genetic 
data, especially given the great intrageneric diversity in 
ecological roles and pathogenicity of bacteria. It, for ex- 
ample, remains unknown why bacteria of the genus Cor- 
ynebacterium are so prevalent in adults, but almost 
absent in chicks. However, some inferences can still be 
made. The seven OTUs shared between chicks and 
adults may be beneficial or commensal and therefore 
retained in the gastrointestinal microbiota as the hosts 
age. For example, two of the shared OTUs were 
Lactobacillales species, many of which are known to com- 
petitively exclude pathogenic bacteria and increase anti- 
body levels, thus increasing immunity to pathogens 
[44,45]. Similarly, another shared OTU belonged to the 
genus Escherichia. Escherichia bacteria are common com- 
mensals in the gastrointestinal tract [46], which are known 
to rapidly colonise the gut of young birds [18,19,24]. 

The data generated in this study will provide new op- 
portunities to investigate the causes and consequences 
of variation in bacterial assemblages in a wild bird spe- 
cies. Although much is known in domesticated bird spe- 
cies, relatively little is known about the acquisition of 
gastrointestinal bacteria assemblages in wild birds. For 



example, assemblages are known to vary with both ex- 
ternal factors, such as nutrition e.g. [47] and environ- 
ment e.g. [48], and host-related traits such as genotype 
e.g. [49], body condition e.g. [48], immune system e.g. 
[50], and sex and mating behaviour e.g. [13]. Our data 
will allow us to build on these studies with large experi- 
mental data sets to explore how specific bacterial species 
are acquired (e.g. via cross-fostering experiments to as- 
certain how environment and genotype affect microbial 
assemblages in chicks) and how they impact on host 
condition and fitness (e.g. whether individuals that host 
different bacterial assemblages cf. enterotypes: [51] differ 
in condition or reproductive success). 

Conclusions 

The striking difference in bacterial assemblages between 
chicks and adults suggests that despite sharing the same 
nesting environment and being fed regurgitated food by 
their parents, others factors affect the acquisition of clo- 
acal bacterial assemblages in kittiwakes. Although sev- 
eral previous studies in both domestic e.g. [17,19,20] and 
wild birds e.g. [10,24,25] have demonstrated age- 
differences in gastrointestinal bacteria, we have done so 
at the community level rather than the single-species 
level. This approach has highlighted how strongly as- 
semblages can differ between age-classes, which was not 
as apparent in previous studies using only culture-based 
methods that target a limited number of bacterial spe- 
cies but see [10]. As gastrointestinal bacteria have im- 
portant functions in digestion, immune functions and 
general health, the characterisation of the acquisition of 
bacteria comprises a crucial component in achieving a 
more comprehensive understanding of the causes and 
consequences of variation in bacterial communities in 
wild animals. 
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